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A physical model for parton densities in hadrons, based on Gaussian momentum 
fluctuations of partons and hadronic baryon-meson fluctuations, is presented. The 
model has previously been shown to describe proton structure function data, and is 
now apphed to sea quark asymmetries and shown to describe the d — u asymmetry 
of the proton. By considering fluctuations involving strange quarks, the model gives 
an asymmetry between the momentum distributions of s and s, which would reduce 
the significance of the NuTeV anomaly. 

1 Introduction 

Asymmetries of sea quark distributions of the nucleon has for quite some time been 
an intriguing problem. For the part of the nucleon sea arising from gluon splittings 
g — > in perturbative QCD, symmetry is expected in the distributions of quarks and 
antiquarks, i.e. q{x) — q{x), and also u{x) — d{x). Conventional parameterizations of 
quark momentum distributions assume these symmetries also for the ^-distributions at 
the start of the perturbative QCD evolution. However, for these sea distributions arising 
from the non-perturbative dynamics of the bound state nucleon experiment has shown 
that there are asymmetries between u and d . There are also no symmetry arguments to 
prevent asymmetries between quark and anti-quark. This is of great interest, especially 
in connection to the NuTeV anomaly, where the value of sin^ 6w was found to differ 
from the Standard Model fitted value by almost three standard deviations 2 . The 
anomaly can however, at least in part, be due to an asymmetry between the momentum 
distribution of s and s in the nucleon sea |31 E] ■ 

Here, we report on recent progress (to be more comprehensively presented in [S]) 
to understand these asymmetries based on a simple and phenomenologically successful 
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Figure 1: Illustration of the processes (a) probing a valence parton in the proton and (b) 
a sea parton in a hadronic fluctuation (letters are four- momenta), (c) shows the resulting 
parton distributions at the starting scale Qq. 

model for the parton distributions in hadrons. In particular, we have shown P" that 
nucleon fluctuations into \A.K), where the s quark is in the heavier A baryon and the s 
is in the lighter K meson, gives a harder momentum distribution for the s than the s. 
This would reduce the NuTeV anomaly to about two standard deviations. 

The model 

Previously, we have presented a physical model giving the momentum distributions of 
partons in the nucleon [H], as illustrated in Fig.^ The model gives the four-momentum 
fc of a single probed valence parton by assuming that, in the nucleon rest frame where 
there is no preferred direction, the shape of the momentum distribution for a parton of 
type i and mass rrii is then taken as a Gaussian 



which may be motivated as a result of the many interactions binding the parton in the 
nucleon. The width of the distribution should be of order hundred MeV from the Heisen- 
berg uncertainty relation applied to the nucleon size, i.e. Gi — l/dN- The momentum 
fraction x of the parton is then defined as the light-cone fraction x — kj^/p+. In order to 
obtain a kinematically allowed final state, we impose the following constraints. The scat- 
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tered parton must be on-shell or having a time- like virtuality, i.e. have a mass-squared in 
the range < < {W is the invariant mass of the hadronic system). Furthermore, 
the hadron remnant r is obtained from energy-momentum conservation and must have 
a time-hke virtuahty. These constraints also ensure that < x < 1. 

Using a Monte Carlo method these parton distributions are integrated numerically 
without approximations. The normalization of the valence distributions is provided by 
the sum rules dx Uv{x) — 2 and dx dv{x) = 1, to get the correct quantum numbers 
of the proton (and similarly for other hadrons) . The gluon normalization is given by the 
momentum sum rule dx xfi{x) = 1, where the sum also includes sea partons. 

To describe the dynamics of the sea partons, we note that the appropriate basis for the 
non-perturbative dynamics of the bound state nucleon is a hadronic quantum mechanical 
basis. Therefore we consider hadronic fluctuations, e.g. for the proton 

\p) = aobo) + ap,rbvr°} + a,„|7i7r+) + . . . + qak |AX+) + . . . (2) 

Probing a parton i in a hadron H of such a fluctuation (Fig. ^p) gives a sea parton 
with light-cone fraction x — xh Xi of the target proton, i. e. the sea distributions are 
obtained from a convolution of the momentum K of the hadron and the momentum 
k of the parton in that hadron. The momentum of the probed hadron is given by a 
similar Gaussian as Eq. but with a separate width parameter gh- The kinematical 
constraints to be applied in this case are mf < p < xhW^ and that the remnants 
(see Fig. n]]) have time-like virtualities. Here xr ~ -M_y/(Mbaryon + -Mmeson), giving a 
harder spectrum for the heavier baryon than the lighter meson, for details see |3|. The 
normalization of the sea distributions is given by the amplitude coefflcients a. These are 
partly given by Clebsch-Gordan coefflcients, but depend primarily on non-perturbative 
dynamics that cannot be calculated from first principles in QCD and are, therefore, taken 
as free parameters. 

This model results in valence and sea parton x-distributions as shown in Fig. ^c). 
These apply at a low scale Qoi ^'^'^ distributions at higher Q'^ are obtained using 
perturbative QCD evolution. The distributions shown in Fig.^c) is the result of the fits 
described below. 

Comparison to data and parameter fitting 

After QCD evolution, the proton structure function F2{x,Q^) can be calculated and 
the model parameters fitted to data from deep inelastic scattering. This results in the 
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Figure 2: The proton structure function F2{x, Q^) from our model compared to HERA 
HI data 7 . 

parameter set 

cr„ = 180 MeV, ad = 150 MeV, Gg = 135 MeV 

GH = 100 MeV, a^^^ = 0.06, = 0.6 GeV^ (3) 

where a^ea fraction of the proton momentum carried by sea quarks at the scale 

Qg. Such inclusive data can only be used to determine the overall normalization aj^^^ 
for the dominating light quark sea from fluctuations with pions in Eq. (j^J. The model 
reproduces F2 data well in view of the its simplicity, see Fig. [5] and [H] . 

The model also reproduces the observed asymmetry between the u and d distributions 
as a result of the suppression of fluctuations with a tt^ relative to those with a 7r+ , since 
the former require a heavier baryon {e.g. A"'"^). With only one additional parameter the 
model gives a nice fit to data on this asymmetry (see Fig.O, the parameter in question 
being the normalization of the |riTr+} -fluctuation relative to the |p7r°)-fluctuation. The 
value of this parameter turns out to be a^^/a^^ = 1/2, which might seem surprising 
in view of the fact that from isospin the relationship should be 2 : 1. However, these 
fluctuations implicitly include the effects of heavier fluctuations like IAtt). The SU(6) 
Clebsch-Gordan coefScient for |A++7r^) is much larger than the coefficient for |A"7r+), 
giving a larger fraction of u relative to d, thus mimicking a larger Iptt*^) fraction (see e.g. 

To fix the normalization of the strange sea in the model, we assume that all fluc- 
tuations including strange quarks (such as \h-K*), \YiK)) can be implicitly included in 
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Figure 3: (a) d[x)/u{x) and (b) xd{x) — xu(x) obtained from the model compared to 
data from the E866 experiment [S]. 

the |Ai4r)-fluctuation and make a fit of the resulting strange sea ((s + s)/2) to data 
from CCFR jTOI. The result is shown in Fig. 01 (left plot), giving a |Ai^}-normalization 
O-AK /olntt ~ 1/5 such that dx{xs{x) + xs{x))/ dx{xu{x) + xd{x)) ~ 0.5, in agree- 
ment with standard parameterizations. 



The NuTeV anomaly and an asymmetric strange sea 

In the NuTeV experiment 2 , the value of sin^ 9w was extracted from neutral and charged 
current cross-sections of neutrinos and anti-neutrinos. The value they find differs by 
about 3(7 from the value obtained in Standard Model fits to data from other experiments: 
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Figure 4: Left: CCFR deep inelastic scattering data ^Hl on the strange sea distribution 
{xs{x) + xs{x))/2 in the nucleon at different compared to our model based on \AK) 
fluctuations (with normalization 's/u'= (s -I- s)/{u + d)) at Qq and evolved to larger 
with perturbative QCD that also adds a symmetric perturbative ss component. Right: 
The strange sea asymmetry s~{x) = xs{x) — xs{x) (at = 20 GeV^) from the model 
and combined with the function F{x) accounting for NuTeV's analysis 
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sin^ 9^'^"'^ = 0.2277 ± 0.0016 while sin^ 6*™ = 0.2227 ± 0.0004. A number of possible 
explanations (31 for this discrepancy have been suggested, both in terms of extensions to 
the Standard Model and in terms of effects within the Standard Model. One explanation 
of the latter kind would be if the momentum distributions for s-quarks in the nucleon 
differs from that of s-quarks. In our model, such an asymmetry arise since the s quark 
is in the heavier A baryon and the s in the lighter K meson, giving a harder momentum 
distribution for the s than the s. 

In Fig. ^ (right plot) we show the resulting asymmetry (x) — xs{x) — xs{x) and its 
combination with a folding function F{x) provided by NuTeV jll| to account for their 
analysis and give the shift in the extracted value of sin^ 0w- We obtain the integrated 
asymmetry 5*" = dxs~{x) = 0.00165, and the shift Asin^^vK = Jq ^~ = 
-0.0017. Thus, the NuTeV value would be shifted to 0.2260 which is only 2.0cr above 
the Standard Model value, leaving no strong hint of physics beyond the Standard Model. 
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